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ABSTRACT 


We have computed the path- i ntegrated gain of parallel pro- 
pagating whistlers driven unstable by an anisotropic distri- 
bution of relativistic electrons in the stable trapping region 
of Jupiter's inner magnetosphere. The requirement that a gain 
of 3 e- foldings of power balance the power lost by imperfect reflec- 
tion along the flux tube sets a stably-trapped flux of elec- 

10 -4 -2 -1 . 

trons J* = 4 X 10 L cm sec which is close to the non-relativistic 
result. Comparison with measurements shows that observed fluxes 
are near the stably-trapped limit, which suggests that whistler 
wave intensities may be high enough to cause significant dif- 
fusion of electrons accounting for the observed reduction of 
phase space densities. A crude estimate of the wave intensity 
necessary to diffuse electrons on a radial diffusion time scale 
yields a magnetic field fluctuation intensity of 
I_ = 1.5 10 /o))^L® ^watts m ^Hz ^ as a lower limit. 
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INTRODUCTION 


RFi^ROnUCIBlun'^ 

ommkal pag:: IS Po<« 

Decreases in observed phase space densities of energetic 
electrons in Jupiter’s inner magnetosphere have led to the sug- 
gestion that pitch angle scattering of the particles by wave 
microturbulence is an active loss mechanism (Fillius et al., 

1975; Mcllwain and FilliuS/ 1975; Baker and Van Allen, 1976). 

A likely suspect is the electron whistler wave or R mode, the 
theory of which when applied to earth’s environment has been 
successful in explaining electron precipitation losses in the 
stable trapping region of the inner magnetosphere (Kennel and 
Petschek, 1966; Lyons et al., 1972). 

Radio observations at decimeter wavelengths of the synchro- 
tron emissions from Jupiter's radiation belts indicated long 
before Pioneer 10 that the electron distribution was highly 
anisotropic and that most of the energetic electrons were con- 
fined to the magnetic equator by a pancaked pitch angle dis- 
tribution (Roberts and Komesaroff, 1965; Thorne, 1965). Pioneers 
lO and 11 confirmed this expectation and also verified the theory 
that the irrmediate source of the radiation belt electrons was 
inward radial diffusion which, conserving the first and second 
particle adiabatic invariants, would flatten the pitch angle 
distribution. Since induced emission of whistlers is a con- 
sequence of anisotropic distribution, the hypothesis that whistler 
turbulence is responsible for observed losses is an attractive 
one both because of the ample growth rates possible and the 
theoretical simplicity of the instability. 


The. first hint that non-synchrotron associated losses were 
present at Jupiter was given by Stansberry and White (197^) 
before Pioneer 10, They set up a radial diffusion model for 
electrons to compute the strip-scan brightness and flux den- 
sity spectrum of synchrotron emissions and found that an ad 
hoc loss process was needed at low L values for a reasonable 
fit with radio observations, Coronlti (197^), in deriving a 
coiTjprehensi ve theoretical model of radiation belt electron fluxes, 
employed the stably-trapped limit concept for whistlers using 
the best pre-Pioneer 10 values for relevant parameters. The 
data from Pioneers 10 and 11 have narrowed the range of param- 
eter space (in particular, the cold plasma density and aniso- 
tropy) available for theoretical models and the observed losses 
invite a reexamination of the relativistic electron-whistler 
interaction. 

The purpose of this paper is to detail some of the con- 
sequences of relativistic electrons in stably-trapped equilibrium 
with parallel-propagating whistler waves. Approximate scaling 
laws for the stably-trapped electron flux and equilibrium wave 
intensity are derived. For simplicity of analysis and clarity 
of content, the major restrictions to our model are the following: 
1) we treat the waves as generated locally and travelling strictly 
parallel to the ambient dipole magnetic field; 2) all reso- 
nant electrons are ultrarelativistic in that the total energy 
is proportional to particle momentum, E==pc. 

In Section II the equatorial growth rate for whistlers 
is derived for a distribution modeled as f(p) “ p” sin^9 . 
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The logarithmic gain for maximally-amplified waves is computed 
by an approximation to the path-integrated growth rate along 
a flux tube and the stably-trapped limit is defined as that 
which will produce 3 e-foldings of power. A similar model for 
the distribution was used by Liemohn (1967) to compute whistler 
amplification by relativistic electrons at earth. Schulz and 
Vampola (1975) considered an artificial radiation belt produced 
by the beta decay of nuclear-fission debris and computed the 
relativistic stably-trapped limit for a distribution with an 
angular distribution similar to ours but with an exponential 
energy dependence. Our procedure and results are consistent 
with the above authors. The stably-trapped limit of relativistic 
electron fluxes is found to compare closely with non-relativistic 
Kennel-Petchek theory. Comparison with measurements shows that 
observed fluxes lie near the limit, which lends support to the 
idea that whistlers are active. In Section III we briefly treat 
the aspects of quasilinear diffusion of relativistic electrons 
and estimate a level of wave intensity that will support dif- 
fusion (losses) on a time scale comparable to radial diffusion. 


S, 
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RELATIVISTIC ELECTRON-WHISTLER INTERACTION 


Equatorial Growth Rate 


If we restrict our attention to whistlers propagating paral- 

A 

lei to the local magnetic field B^z, the dispersion relation 

for R modes with » 0 is (Lerche, 1969) 


K(k. «)) « 0 


kI ^+1 i V^TTe^ r ^ 

2^ Zi moo i ^uoY - kp,, 


uoY " kPw + n 


* l^Px 


J^fp Iffi] 

muoYl^-*- bpjJ 


The ^ represents a sum over particle species (species labels 

on all quantities are implicit), y = T 1 + P^T^^c^ p , m is the 

qB ^ 

particle rest mass, and 0 “ particle distribution 

function fQ(p) is composed of a cold plasma background and a 
small relativisxlc component f, then defining t = p/mc and 

P . B 

X « cos 0 = — 5 — , equation (1) can be written in spherical 
P°o 

coordinates, assuming azimuthal symmetry in p-space, as 


l2 2 

k c 


= 1 + 


.V 


i»l dx r V « w 


•J ^ , /U)V +nMbt bxVfOY tJj ' ' 


where 


2 

“*p ” m 


and n. 
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OKKUSAi- 

for the electrons. If the wave number has a small imaginary 

- K. 

part, k * k^+ikj, we can approximate kj integrating 

around the indented pole of the integrand for the electron con- 
tribution yields 

° X • — 7 

‘"Y • flee 

The integration is taken along the resonance contour x ■ — — 

and t^ represents the minimum momentum (magnitude) satisfying 
wYo " flee 

-1 ■ — kc T ^ * subscript on the real part of the wavenumber 


as been dropped rkj.-»kJ and 


k^c^ , 

— s— ■ n ■ I +• 


.^cl 


The resonance curve in momentum space is a hyperbola, plotted 
in Fig. 1, satisfying the relation 



The resonance curve crosses the axis at the point tj^ = 
and asymptotes to the angle cp ■ sin”^ The passage to the 
non-relativist ic regime requires n»l, t«l and the resonance 
curve flattens to approximately a straight line. Equation (6) 
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V. V« 1 . H -o 

then gives the familiar result -t^ » tg where — " n *“ n 

for 0 ) « 1 and Eg * (Kennel and Petchek, 1966). Relativistic 

particles can interact with whistlers either by Ooppler-shI ft Ing 

the wave frequency up to the local electron gyrofrequency or 

by lowering their gyrofrequency sufficiently by a mass Increase, 

thus the resonance curve departs from a straight line. Induced 

emission or absorption of the whistler depends on the local 

derivatives of the distribution function along the resonance 

curve. The growth or damping of a mode depends on the net energy 

contribution to a wave from the particles,' For a distribution 
M 

such as f«sin 9, a given frequency may be unstable by equation 

(3) but for those particles with p_ ■ 0, ^ *0 and con- 

^ 0-n/2 

sequently those particles are energized by the growing mode. 

For application to energetic electrons in Jupiter's magneto- 
sphere it is convenient to use the ul trarel at 1 v i,st i c approxi- 
mation where t»l. In that case the resonance curve can be 
approximated by x^^[l-n/t] valid so long as t^ V? n ^ * 

If the relativistic electron distribution obeys a power law 

In momenta (energy) we can model the distribution at the mag- 

_ 

netic equator as f(t) ■ where M and N are the pitch angle 

and spectral indices of the distribution and B is a normali- 
zation. A convenient normalization is in terms of the omni- 
directional integral flux. For relativistic particles where 
Y ■ ,/l + t^ ^ t we have 




r 2 

such that J(>y) ■ j dvj cf(t). Since the growth rate depends 

Y 

on the distribution above t^ by (3), If we consider frequencies 
such that t^ lies In the domain where (7) Is valid, we do not 
have to specify the behavior of the distribution at low energies. 
When (7) Is Inserted Into the growth rate and the change of 
variable x ■ ^1 -Q/t] Is made we have the following expres- 
sion for the logarithmic gain scaled to a Jovian radius at the 

« 

magnetic equator: 




,N-1 


■1 F(N, M, u) 


( 8 ) 


where 


/M + 3N 


F(N, M, u) . (N- dx(u-x) 


N-1 




pTniT 




M 


* “X^)^[u(N + 2)f 1 -x^) +Mx(l - ux)] (9) 


and u - ^ is the normalized phase velocity of the wave. The 
growth rate In the relativistic approximation is dependent more 
fundamentally on the refractive index (actually phase veloc- 
ity) rather than frequency. 

If the quantity tp * n/n is considered as an effective 2 

^ TT eRj 

resonant momentum, we can Interpret equation (8) such that 
is the electromagnetic coupling of the wave to the distribution, 
is the number "density" of particles with momenta 
titg, and F(N, M, u) is the kinematical resonance integral 
over the curvature of the distribution. Figure 2 is a plot 
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of F(N, M, u) for N ■ 3 and M ■ 0, 2, 4, 6, 8. Negative values 

of correspond to growth. At low frequencies, «)-*0, the 

N 

growth rate goes to zero falling off like tjJ . At higher fre- 
quencies where u-0,3t the growth rate turns over and goes to 
zero since at a high enough phase velocity as much energy Is 
taken from the wave as Is given, the resonance curve In Fig. 

1 becoming Increasingly concave. 

If the cold plasma density Is known the growth rate can 
be computed as a function of frequency through the transforma- 
tion — - n(«)), A basic role of the cold plasma density Is to 
determine a range of possible values for the refractive Index: 

"mAX^'’’ ■ ^ ^ ”MIn("p 2^ce)» '^MIN o^^^^rrlng at the lesser 

of (ttp O'* range specifies a phase velocity window 

which restricts the values of u In Fig, 2. For Jupiter, Frank 
et al. (1976) have reported cold proton densities, which have 
been used to plot Fig, 3. The cross-hatched region represents 
the accessible values of phase velocity at each L = on 

the magnetic equator. By comparison with Fig. 2, using repre- 
sentative values of M » 4 and N • 3 for the pitch angle and 
spectral indices, there exists a range of unstable frequencies 
for 1 ^ L < 1 2, 

In the non-relatlvlstic regime one can determine the mar- 
ginally stable frequency and the corresponding resonant par- 
ticle energy above which whistlers are unstable. For relativistic 
particles, since the resonance curve is hyperbolic, the mar- 
ginal stability point Is model -dependent and a function of M 
and N, We can estimate ^ as follows. Choosing n « 3 and 
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V i'V * ii 


M ■ 4, the marginally stable refractive ri|^ 5 ••2, 

E|^ ^ Is non-relat ivistic and all relativistic electrons resQ; 
nate with unstable whistlers. If n^ ., » 0(1), t^ - t^ c t ^ i « 

2 1' ■ O n# ^ a fli^ C ^ • 

/ \ n.i. 

Taking «» « 1 gives tj^ g » S 'V 


w. 


^ ("M.S.- ’) “T!^ 
U)«- o 


pc 


essentially the non-relativistlc result. 


b) Path Integrated Growth 


( 10 ) 


The actual logarithmic power gain of a ducted whistler 

fS? 

requires the path-integrated quantity 6 *-2j k|(s)ds along 
a flux tube. With certain approximations the integration can 
be done analytically and the result is good at least to the 
accuracy that fluxes are measured. We treat the cold plasma 
density as constant along the flux tube and require equal con- 
tributions to the gain from above and below the magnetic equator. 
As the wave convects along the tube from one hemisphere to the 
other, the phase velocity changes as demonstrated in Fig. 2 
by either of the arrows A or B depending on the phase veloc- 
ity (frequency) at the initial point. Liemohn (1967) found 
that waves which were locally damped at the magnetic equator 
could have a net positive gain because of greater growth con- 
tributions at higher latitudes. This situation is demonstrated 
by the arrow A for the curve M » 4. However, the stably-trapped 
flux is determined by the waves with maximal gain, which case 
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la represented by the arrow B; the situation marries the maxi- 
mum of F(M, N, u) with the maximum flux of resonant particles 
at the equator. But from inspection of Fig. 2 , if we choose 
N ■ 3 and M ■ 4 the excursion due to F(M, N, u) is small pro- 
vided we consider phase velocities such that u(u)) < .35. We 
thus approximate F(M# N, u) as constant in this frequency range 
and having the constant value F(N « 3, M ■ 4, u = 0) « -5/16 
(dotted lines) along the flux tube. The path integration then 
simply involves the scaling of the other factors in equation (8) 
along the flux txibe. The coupling scales like off the 

equator. For a sin^Q distribution the flux scales like (B^/B)^^'^^ 
(RoeJirtrer, 1970). If n “ 0(1), tj^ = n/n scales like (B/B^) . 

The growth rate along the tube is then approximately 

B ilt2N B 5 

ki(s) - (-^) 2 k^(EQ) =(-^) k^(EQ) (11) 

s 

Then if G ■ - 4/^”'“^k^ (s) ds and using the harmonic approximation 

2 2 2 2 2 2 
B/Bq k 1 + s /s^ for a dipole magnetic field where s^ » Rj 

we have 

y - (M^-2N) 

G -• -4k^(EQ)LRj^/^"'^^dy(l + y^) ^ (12) 

When y = 1 we are at the limit of the harmonic approximation, 
■*max 

but since the integrand has decreased significantly, we extend 
the upper limit to infinity and for N = 3, M = 4 the integra- 
tion yields G = - 4k^RjL(0.2) . Thus 
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G,h=3, M=4)= 


r- -i2 
Y 

fl/n 


(13) 


Q 

We note that in scaling t_ - (=— ) we required n = 0(1). If the 

K O 

^ B 3 /2 

opposite extreme prevailed, n>> 1, then tj^'^' (g— ) ' and 
B 6 o 

(s) - (-^) k. (EQ) , with the integration providing G = - 4k.R-.L(0 

X o X X ij 

a slight difference. All quantities in equation (13) refer 
to the equator and (13) is valid for frequencies such that 
u(ai) < 0.35. 

Generalizing this procedure for arbitrary values of M and 
N and approximating F(N, M, u) ==F(N, M, u = 0) we have 




r(^4^) r(- 


^(N- 1) 


lA) r/ M-i-2N- 1 > 2 

2 ^ 2 ’ J(>Y) 






Y 

H/n 


We note that over the parameter range M=2, 4, N=2, 3, 4, 

G is a relatively weak function of the indices, viz. 

G(N = 2, M = 4) ='(1.27) G(N = 4, M = 2) at the extreme. This 
is due to the compensating effects of, say, larger growth rate 
(M increases) with smal?.er effective path length from equation 
(11). Therefore, over this parameter range we will take equation 
(13) as valid generally. 


18) , 


N-1 

(14) 
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* c. Stabi v-Trapped Limit 


Application of equation ( 13 ) to the concept of the stably- 
trapped limit (Kennel and Petchek, 1966) requires a knowledge 
of the power reflection coefficient, R, for the whistlers in 
the flux tube. If we treat the problem as strictly one-dimensional 
and define the volume emissivity r\ for paral lei -propagating whistlers, 
the equation of radiative transfer ^ *-2k; (s )J_ + <q(s ) yields 
for the intensity 




dsTi(s) exp.-2j^'^^^ds'k.(s')l 
‘^MAX L s J 


1 - Re' 


(15) 


Q 

It is clear that as Re -* 1 , sufficiently high wave intensities 
w.i 1 1 result that can relax the distribution on a time scale 
comparable to that of the particle source. No in situ wave 
measurements have- been conducted at Jupiter, but if we assume 
that as In the Kennel -Petchek theory 5% pf the wave energy is 
reflected, a gain of G = In (l/R) = 3 is sufficient to main- 
tain the stably trapped equilibrium. Equation ( 13 ) then gives 
a stably- trapped flux limit of 



4 ^ , o' O l-U electrons 
cm sec 


(16) 


The non-relativist ic Kennel - Petschek result scaled to Jupiter 
is 


■'i5.R(>ER^ ' ^xio’VV^) = S.Sxlo'V'* (17) 


t-iittpIIIlY of ftiD 

reproouc&h^ 

OUIU’.NAI. 
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which does not differ greatly from the relativistic result. 

This fact is borne out by inspection of (8). The relativistic 

aspect is manifest only in the resonance integral F(N, M, u) 

and If the distribution is not too different topologically in 

the non-relativi St ic regime, F(N, M, u) should not vary a great 

deal. This can be further motivated by assuming a power law 

distribution In non-relativistic energy. If we take the dis- 

• !■ 

tributlon with = T as f(v) » where M and N are the 

pitch angle and spectral Indices, a npn-relativistlc analysis 
yields the result equivalent to equation (8) for tu « 1 of 

2_o f T I T NN"1 
TT 

•rj “ " 


k.R. - . V/— R . - _. Fn.r(N, M. «) 


B. 


(18) 


where 


and 


, , 1 w2 1 /'0c>^ 


* -M(N- 1)^ 


rl 

^M + 3\ 


V 2 ) 

r( 

+ 2\ 


< 2 J 


M(1 


- iu)_ 


c- M+K 2N-1„ 
. J d0 s I n 9 cos 9 


09) 


1 f uj « 1 and we neglect the second term in (19), then for N « 3 
and M = 4, * “7^ compared to F =*■ ■ taken previously. 

Equation (16) describes a limit which should not be exceeded 


statistically by a flux of electrons at .any energy. If the 
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flux J(>y) is known and obeys a power law, the distribution 
should depart from the power law and harden considerably at 
energies below the transition momentum given by 


t 


R 




1 


r “1 

TTHT 


J(>v) . 


Y 




( 20 ) 


Since further extrapolation of the distribution by a power law 
to lower momenta would produce excessive amplification of the 
whistlers. 

In Fig. 4 we have plotted J*(>tj^) for comparison with observa- 
tion. For L£l2 measured electron spectra have spectral indices 
from 3-3.5 at high energies and pitch angle indices from 2-4 
(Van Allen et al., 1974; Baker and Van Allen, 1976). We have 
taken the representative values of N = 3 and M = 4 throughout. 

The circles are equatorial electron fluxes of 5 Mev electrons 
from an empirical formula given by Mcllwain and Fillius (1975) ; 
the diamonds are 5 Mev fluxes reported by Baker and Van Allen 
(1976) with error bars (D.N. Baker, private communication) . 

The theoretical uncertainty is more sensitive to the basic 
model we have chosen rather than to M and N values. We have 
taken the most efficient situation of parallel-propagation for 
which the growth rate is maximum (Kennel and Petschek, 1966) ; 
the finite kj^ effects should raise J*. A reflection coefficient 
of 0.5% rather than 5% raises J* by a factor of 1.8. Another 

f 

possibility is that centrifugal forces at large L confine the 
cold plasma to within ±1 Rj of the equator as suggested by loannidis 
and Brice (1971) . The effective path length is then just 2 Rj 


and for N « 3, M « 4 the modification to the limiting flux becomes 
9 -3 -2 -1 

j* '\» 8 X 10 L cm sec which is also plotted in Fig. 4. Within 
the experimental emd theoretical uncertainties involved, the 
data lie near the stably-trapped limit and suggest that whistlers 
may be active and that the transition energy is approximately 
5 Mev at low L values. 
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III DIFFUSION THEORY 

t 

a. Homogeneous Quasi -Linear Theory 


If we consider a homogeneous plasma with paral lei -propagat i ng 
whistlers described by a one-dImens ional electric field spectral 
density 6(uj) for waves travelling in just one direction, then 




^'TOT “ 


du)^!(uj)j 


the factor of two includes the waves travelling in the oppo- 
site sense. The quasi linear diffusion equation for the reso- 
nant electrons can be written 


where 


bT / 


Q.L. 


J 

bt 


( 21 ) 


I I'^gl IV 




m c 


nt 


g 

nt 


z 


l^bt 


jlL+h t - 1 


/ '■i 


V 'z 


( 22 ) 


and tj^, t are unit vectors. Since a relativistic particle 


interacts with two waves going in opposite directions, the sum- 

^ce 

nation is over the frequencies satisfying ' 


and 


k I c 


j'ections, the 
' =■ 


for the R mode travelling in the + direction 


along The combination in the first bracket is proportiona' 


r, I'I'Y < ‘.\j' itllij 

. i»: I . = . ^ ^ 
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to the incremental growth rate (see equation (1)) and has a 

sign, 0 ■ +1 (-1), corresponding to the particles giving (taking) 

energy from the wave. The diffusion current is thus proportional 

to the spectral density and incremental growth rate and, if 

we consider the interaction with just an upward travel 1 i ng wave, 

has the direction given by the angle K (see Fig, 5), where 
J ^ 

♦+ " * c - nv eos^ e ' '' 




(23) 


in general, is directed such that those particles giving 
energy increase their pitch angle, 9-n, and those receiving 
energy decrease their pitch angle, 9-0, If n cos 9«1, 
t^“*tan - tan 9]+y(1+o), which describes pitch angle dif- 
fusion along approximately iso-energy surfaces, since 
t^ + t^ = constant leads to = - tan 9. At 9 » tt/2 = tan’ n, 
but these particles react equally to downgoing waves and the 

A 

net current is in the tj^^ direction. Thus 90° pitch angle par- 
ticles should random-walk along the t^ axis to higher energies. 

The features of the diffusion are sketched qualitatively in 
Fig, 5. The locus of points where the incremental growth is 
zero defines a cone inside of which the diffusion is generally 
directed outwards. We note that non- rel at i vi st i cal 1 y , i f we 
use the resonance condition 1 cos 0 * 




tan’ (u) - 1 )tan 9]+^(l+o) 


(24) 
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and for low frequencies, ua « 1 , pitch angle diffusion results, 
b. Radial Diffusion Particle Source 

The Immediate source of relativistic electrons In the Inner 
Jovian magnetosphere has been Identified as Inward radial diffusion 
resulting from third adiabatic Invarlation violation (see, e.g,, 
Simpson et al., 1974). As a source term the process can be written 



where the derivatives, are taken at constant first and 

second Invariant, M and J. 1 f we neglect synchrotron losses 
and bounce average, <), equation (21), the evolution of the 
distribution under combined radial diffusion and quasi linear 
momentum- space diffusion by parallel-propagating whistlers is 
described by 

^ s (26) 

bT bryR.D. 

If electric field intensities are small the first term domi- 
nates, but If the distribution is sufficiently unstable to whistlers 
and in a stably-trapped equilibrium, quasi linear diffusion must 
occur on a time scale comparable to that of radial diffusion. 

The solution of (26) In the steady state Is formidable, but 
we can extract a crude approximation of the bounce-averaged 
<C(uj)> in the stabl y- trapped regime. The radial diffusion coefficient 
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can be modeled as °LL ■ “o'-” where — lo’ ^^sec’ ^ and a -4 (Barbosa 
and Coronltl, 1975 ). Thus we estimate Writing 

Dq l 

(21) In spherical coordinates, we estimate the Q.L. relaxation 

2^2^2 |V^-v_ 

time as 2 2 <Vft(m)) * crudeness of the approxi- 

mation Is manifest In the neglect of the pitch angle dependence 
of tq l also by the neglect of the derivatives In (22) 

and ( 25 ) which to some degree are compensatory. Equating charac- 
teristic times yields 

<Vge(u,)> - 

If Iv - V 1 c amd t = t_ = iT/n, then letting 2ttv « oj we have 

I 9 21 R 

for the magnetic field intensity of modes travelling parallel 
to the field in one direction 

♦ 

<I„(v)> = 27 T<n^I (w) > 'V 1.5 X W^TS ^28) 

m Hz 

<I > is a lower limit such that intensities much lower than 
B 

(28) will not produce significant decreases in phase space den- 
sities evolving under pure radial diffusion. 

We Ccin estimate typical fluctuation field strengths from 

n 2 

(28). If the bandwidth Av 'v. v 'v then ^Av< 3 (v)> 

and = .5 myL”** for a = 4. Such fields should 

be detectable by future spacecraft to Jupiter. 
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tv DISCUSSION 

The evaluation of the limiting flux assumed whistlers were 
generated locally and did not propagate across L-shells. With- 
out more detailed knowledge of the cold plasma distribution, 
ray path computations would be speculative. However, the fact 
that observed fluxes are lower than the limiting flux inside 
of L<4, where phase space losses are still apparent, suggests 
that whistlers may propagate inward. 
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curves plotted In normalized momentun space 
- The solid line is the resonance with 
a whistler (Q ■ 10, n * 1.94) travelling along the 
field; the dashed line with a whistler (n ■ 4, n ■ 1,53) 
travelling anti-parallel. 

Plot of the resonance function F(M, N, u) for M - 0, 

2, 4, 6, 8 and N 3. Growth occurs for negative values 
of F(M, N, u) increasing with M. The arrows repre- 
sent convective changes of F along the flux tube. 

Accessible phase velocities for whistlers using Frank 
et al. (1976) cold plasma observations. 

Plot of the limiting flux J*(>tj^). The circles are 
5 Mev equatorial fluxes of >^,cllwain and Fillius (1975) ; 
the diamonds are 5 Mev fluxes of Baker and Van Allen 
(1976) with error bars. The dotted line is the modi- 
fication to J* from centrifugal effects. 

Qualitative view of diffusion in normalized momen- 
tum space arising from interaction with whistler pro- 
pagating parallel to only. The small arrows give 
only the direction of the diffusion. The cone edge 
is defined by the integrand of (3) being zero for N ■ 3 
and M = 4. 
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